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Abstract 

A statistical model for decay and formation of heavy hadronic resonances is formulated. The 
resonance properties become increasingly uncertain with increasing resonance mass. Drawing on 
analogy with the situation in low-energy nuclear physics, we employ the Weisskopf approach to the 
resonance processes. In the large-mass limit, the density of resonance states in mass is governed 
by a universal Hagedorn-like temperature Tjj. As resonances decay, progressively more and more 
numerous lighter states get populated. For Tjj ~ 170 MeV, the model describes data for the 
hadron yield ratios at the RHIC and SPS energies under the extreme assumption of a single heavy 
resonance giving rise to measured yields. 

PACS numbers: 12.38.Mh, 24.85. +p, 25.75.-q 
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Studies of ultrarelativistic nuclear reactions aim at learning on properties of highly-excited 
strongly-interactingmatter and, in particular, on the predicted transition to quark-gluon 
plasma (QGP) jllO]. Yields from those reactions can be described in terms of grand- 
canonical thermodynamic models DDDQ , at temperatures close to those anticipated 



ermoayna 
sition [7[. 



for the transition |7[. Similar temperatures are utilized in the microcanonical model for 
electron-positron annihilation into hadrons 8] . Large fractions of final light-particles within 
those thermodynamic equilibrium models result from secondary decays of heavy resonances 
whose features are, generally, less and less known the heavier the resonances. The lack of 
aiowlcdgc has forced the use of cut-offs on the primary resonance mass in equilibrium models 
Mft analog cn t - offs have been e mploy ed in „ m „dels HflQ. In this 
paper, we explore the possibility of an apparent equilibrium in the reactions stemming from 
sequential binary decays of heavy resonances, with the characteristic temperature describing 
the resonance mass spectrum, rather than being proposed ad hoc. 

The situation of deteriorating knowledge of resonance properties, with increasing reso- 
nance energy and resonance density in energy, reminisces that of resonances in low-energy 
nuclear physics. There, statistical descriptions, in terms of the Weisskopf compound-nucleus 
model and the Hauser-Feschbach theory, have been highly successful 12|. Underlying the 
statistical descriptions is the density of resonance states in energy p(m, q), where q represents 
discrete quantum numbers of the resonances. For the hadronic states, a universal tempera- 
ture in terms of the density of states, = dlogp/dm, has been originally considered by 
Hagedorn Q], in part due to an evidence for a rapid, nearly exponential, increase in the 
number of resonances with energy; see also 3]. To reach a conclusion on the basis of spec- 
trum, however, the counting of resonances should be carried out at constant values of the 
quantum numbers q, particularly in the region of opening thresholds for different q. This has 



not been considered in Refs. 



13 



14j . Nonetheless, in the limit of large m at fixed q/m, the 



universal temperature Th, independent of m, but possibly dependent on q/m, may be ex- 
pected for the resonances due to the lack of any scales that could govern the T^-dependence 
on m. While our model for p will be fairly schematic, similar to the models employed in 
the literature 13, 2], there are generally important questions regarding strong interactions 
that can be suitably asked in terms of p, concerning e.g., besides the T H (q/m) dependence, 
the emergence of a surface tension in the thermodynamic limit. In the context of the phase 
transition, the temperature Tjj may be considered as the temperature for a metastable equi- 
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librium of quark-gluon drops with vacuum, and, as such, slightly lower temperature than 
the critical T c . 

For resonances described by the continuum density of states p, we consider the processes 
of binary breakup and inverse fusion, constrained by detailed balance, geometry and by the 
conservation of baryon number B, strangeness S, isospin I and of isospin projection I z . 
From the two versions of our model, with and without a strict conservation of angular mo- 
mentum (J, J z ), we discuss the simpler Weisskopf version, in this first model presentation. 
The lighter particles (comprising of 55 baryonic and 34 mesonic states) are treated explic- 
itly as discrete states in the model. Our model allows to explore various aspects of the 
system evolution in ultrarelativistic collisions, including formation of hadrons with extreme 
strangeness and isospin, as well as chemical and kinetic freeze-out. 

We first discuss details of the density of states. As a threshold mass, separating the 
discrete states from those in the statistical continuum, we take mj: h = 2 GeV. We assume 
that at low excitation energies the density of states is similar at different q for the same 
excitation energy above the spectrum bottom m g (q). Following quark considerations, we 
adopt 

m g (q) — ciq max (| 3B + S\, 2 I) + as \ S\ , (1) 

where = 0.387 GeV and a s = 0.459 GeV. The coefficient magnitudes have been ad- 
justed requiring that the mass m g (q) from Eq. exceeds the lowest known masses for 
different q; regarding the practicality of sequential decays, we prefer to overestimate rather 
than to underestimate the reference position of spectrum bottom for the continuum, and 
to overestimate its rise with q, to preclude the emergence of any unphysical stable states 
towards the continuum bottom. At high masses m the influence of the ground state mass 
on p should decrease. We eventually arrive at the following density of states employed in 
our calculations: 

p(m, q ) = A ^ [{m - m ' f{m - m J }/T *L (2) 
[{m — m g f(m — m g )} + m%\ 

The role of the factor / is to suppress the effect of m g at high m and we use 

f(m - mg) = —j- — n . (3) 

1 + [(m-m g )/m c th \ 

with n = l. The prefactor of the exponential in Eq. (|2"|). with m r = 0.5 GeV, acts to enhance 
asymmetric (rather than symmetric) breakups for continuum hadrons, leaving the issue of 
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surface tension in the thermodynamic limit open. The normalizing factor A and the power 
a are adjusted, for different assumed values of Hagedorn temperature Th, by comparing the 



low-m cumulant spectra from measurements and from the continuum representation 

N exp = ^(2J+l)6(m-m 4 ), 

i 

N the = J2 / dmV(TO',g). (4) 



<? 

For Th = 170 MeV, as an example, we obtain the prefactor power of a = 2.82. The 
normalization factor A drops out from probabilities for the most common processes involving 
continuum hadrons, with one continuum and one discrete hadron either in the initial or final 
state. 

We now turn to cross sections and decay rates. The cross section for the formation of 
a resonance q in the interaction of hadrons 1 and 2, can be, on one hand, represented as jl^l ]: 

/ \ 1 m i m 2 /"dp f de , . , , l2 

V12 ei(pi) e 2 (p 2 ) y (2tt)' 3 J 2vr 
x (2tt) 3 5 (px + p 2 - p) 2tt 5 (ei(pi) + e 2 (p 2 ) - e) 

■p{m,q) \M qi +q 2 ^ q \ 2 . (5) 



mp*(mi, m 2 ) 

Here, t>i 2 is the relative velocity, e's are the single-particle energies, and \ Ai\ 2 is the matrix 
element squared for the fusion, which is averaged over initial and final spin directions. The 
factor of (2 J + 1), associated with the last averaging, has been absorbed into p. The cm. 
momentum in Eq. (jSJ) is 

p*{m a , m b ) = — [(m 2 - (m a + m b ) 2 ) (m 2 - (m a - m b ) 2 )] 1/2 . (6) 

For the final state in continuum, following geometric considerations, the cross section for 
fusion, on the other hand, is 

<r(qi + Q2^q) = (h I Z1 1 2 I Z2 \\II z )nR 2 , (7) 

where (-||-) represents the isospin Clebsch-Gordan coefficient. The cross section radius is 
taken as that of the fused entity, R ~ tq (m/nid) 1 ^ 3 , with r$ = 1 fm as a characteristic radius 
for a rricL = 1 GeV hadron. Equations (j3J) and (J7|) allow to extract the square of the transition 
matrix element needed for computation of the partial resonance-width. 
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The partial width for decay into q\ and q 2 can be generally represented as 

T(q - ft + fc) = [ -% [ dm[ ft) / dm' 2 ™ 2 ^ ^ 



(2k)* J 1 ei (p) J ' e 2 (p) 
x\M n+q2 \ 2 2n 5 ( ei (p) + e 2 (p) -m) . (8) 

As before, the factors of (2Jj + 1) are absorbed into pj. A resonance in the continuum 
can undergo three types of binary decay: where (i) both daughters are particles with well 
established properties within the discrete spectrum below mjr h = 2 GeV, (ii) one of the 
daughters belongs to the discrete spectrum and the other to the continuum and, finally, where 
(iii) both daughter resonances belong to the continuum. In the case (i), the state densities 
are Pj{m') = (2 J,- + 1) 5{m'- — raj). With the detailed balance relation, \A4 q ^ qi+q2 \ 2 = 



3 
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l-A4gi+g a -«jl , we then get 



mi 777.2 p* (mi, 



r«( ? -g 1 + ?2 ) = """^^'^ (2Ji + l) (2J 2 + 1) |A4^ J 2 

7rm 

_ lT T j j (2J 1 + l)(2J 2 + l)^ 2 (m 1 ,m 2 )fi 2 

- (/1 J Z1 7 2 J 22 1 1 J J,) 2vr p(m,g) ' (9) 

In analyzing the case (ii), let the particle characterized by q\ belong to the discrete 
spectrum and that characterized by q 2 to the continuum spectrum. From (JHJ), we then find 

?2 



r ( %-g 1 + g 2 ) = (hI zl hU\II z ) W + 1 ) mjy 

Z7T p(m, gj 



x r* (mi ' me2) dp p 3 p(y/m* + m\ - 2m ei ,g 2 ) 
Jo ei ^m 2 + m 2 — 2m ei 

where m C2 = max(m£ h , m s (g 2 )). For ei small compared to m, it is of advantage to represent 
the subintegral density of states as an exponential of the density logarithm and to expand 
the logarithm in (ei— mi). The integration over relative momentum can be thereafter carried 
out explicitly, yielding 

7rm 2 p(m,qj 

Here, m 2 = m — mi and the temperature is defined as 

1 m d logp(m 2 , q 2 ) 
T 2 m 2 <9m 2 

m J 1 | n(m t c h ) n [m 2 - m^fe)]""' ■ ^ 



m 2 T^ ) (mt h ) n + [m 2 - m g (q 2 )] n 
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For (mi,qi) small compared to (m,q), the obvious further possibility is the expansion of 
logp(m2,q , 2) in (mi, gi), with an emergence of the chemical potentials conjugate to q 2 . 

In the case (iii), of both daughters in the continuum with large masses m ; > m c th = 2 GeV, 
the nonrelativistic limit in Eq. (jHJ) is justified. On employing Eqs. (0) and (J7J), the result for 
the partial width is 

r (^(q^ qi + q 2 ) = (1,1 IzIza \\ IIz ) / deep(m-e, qi ,q 2 ) 

Sir p(m,q) Jo 

= /,/,l|//.) '"f l57K 7 g "f ) , (13) 

8n p(m, q) 

where 

™-2m C2 

p(m,qi,q 2 )= / dup ((m + u)/2, q x ) p ((m - u)/2, q 2 ) (l-w 2 /m 2 ) . (14) 

J 2m Cl — m 

To obtain the last expression in Eq. ([13)1 . we have expanded the logarithm of subintegral 
density, with the temperature representing 1/Xi 2 = d[logp(m,q 1 ,q 2 )}/dm. 
The total decay width of an (m, q) resonance is finally 

r(m, q) = J2 rW (? ^ 9i + ? 2 ) + X] rW (9 ^ 9i + 92) + ^ r( is )( g 9l + g 2 ) . (15) 

(?i 92 91 92 91 92 

A moving resonance will live an average time of (r) = 7/T, where 7 is the resonance Lorentz 
factor. 

The formulas above provide the basis for our Monte-Carlo simulations of the resonance 
decay sequences in heavy ion collisions. A resonance follows an exponential decay law 
corresponding to (r). The product properties are selected according to the decay branching 
ratios. Since the parent angular momentum is not tracked in the Weisskopf approach, the 
angular distribution of products is taken as isotropic. During the evolution two resonances 
can fuse with each other, according to the cross section of Eq. (JJJ), if the final state is 
in continuum. If, on the other hand, the state is discrete, the cross section acquires the 
standard form [l5j|, from Eqs. (0) and (jHJ), 

/ , \ 2tt 2 p(m 1 q)T(q -> q x + q 2 ) 

(2J 1 + 1 )(2J 2 + l)^(m„m 2 ) - (16) 

where the width for the spectral function is considered: 

(2J + l)2m 9 2 r ^ (2J+1)T 
P[7U ' q) 7T [(m 2 - + m 2 T 2 ] ~ 2vr [(m - m,) 2 + T 2 /4] ' 1 ' 
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Besides the decay and fusion processes, related by detailed balance, a provisional constant 
cross section a e \ = 5 mb has been assumed for all collisions. 

In our model, we simulate, in particular, the features of the final state of central Au+Au 
collisions at y'i = 130^4 GeV. Following the presumption that the decay sequences will tend 
to erase fine details of the initial state, we push the characteristics of the initial state to an 
extreme, allowing for a single heavy resonance to populate a given rapidity region. In the 
end, when examining the transverse momentum spectra, we find that the resonance decay 
and reformation alone generates insufficient transverse collective energy, indicating that the 
early resonances need to be affected by a collective motion generated prior to the resonance 
stage. This finding is consistent with those in other works 0, Q|- 

Within the single resonance scenario, the local final state reflects the initial quantum 
numbers of a resonance characterized by (m ,g ), where q = (B , S , 1 , I z0 ). After the 
value of the Hagedorn temperature T H is set, the relative yields of particles in the final state 
are, in practice, sensitive only to the relative values of the quantum numbers of the initial 
resonance. We normally impose strangeness neutrality, so that the starting value is So = 0. 
The magnitude of Bq, for a given mo, can be adjusted by using the final-state antiproton- 
to-proton or antiproton-to-pion ratios. The starting isospin, for a given B , can be adjusted 
by using the isospin of original nuclei. However, we find that the initial isospin has only 
a marginal impact on the isospin of individual final particles. This may be attributed to the 
cumulative effect of isospin fluctuations when many particles, compared to /, are produced. 
For specific initial (mo, qo) values, we repeat numerous times the Monte-Carlo simulations 
of the decay chain and recombination, and the results presented here are an average of 
about 10 4 generated event sequences. 

Figure ^ illustrates the average features of an exemplary local system that starts out as a 
resonance characterized by m = 100 GeV and B = 3. The left panel shows the ratio of the 
average maximal resonance mass m max to the initial mass mo, as a function of time, as well 
as the mass asymmetry, the mass difference between the heaviest resonance and the next 
heaviest, divided by the sum of their masses, a<i = (m max — m2)/(m max + m 2 ). Persistence of 
the large asymmetry with time indicates that the heavy resonance decays primarily through 
light-hadron emission. 

In addition, Fig. ^ shows, as a function of time, the abundances of particles (left panel) 
and antiparticles (right panel). The antibaryon abundances freeze out noticeably earlier 
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FIG. 1: Time evolution of average characteristics for a system starting at mo = 100 GeV, Bq = 3 
and Sq = 0, under the assumption of Tjj = 170 MeV. Displayed in the left panel is the ratio of the 
maximal resonance mass to the initial resonance mass, m max /?Tio, as well as the mass asymmetry, 
a>2 = (fftmax — m2)/(m max + ^-2), where is the mass of the second heaviest resonance. These 
results are shown both for the standard system evolution, then represented by the monotonically 
dropping solid and dashed lines, and for the evolution with suppressed fusion reactions, then 
represented by the monotonically dropping dash-dotted and dotted lines, respectively. In addition, 
the left and right panel display the evolution of particle and antiparticle abundances, respectively, 
for the standard evolution, in terms of rising lines. 

than the baryon abundances, in spite of a low initial baryon number relative to the initial 
mass. The higher the strangeness, the later the abundance saturates. This is likely due 
to the fact that the effects of strangeness fluctuation need to accumulate with time; the 
situation would change if we assumed strangeness fluctuations right for the initial resonance 
conditions. 

The calculations have been repeated while suppressing the back resonance fusion reac- 
tions. As expected, in this case the maximal mass and asymmetry decrease faster with time, 
see the left panel of Fig. ^ The abundances (not shown) grow faster and saturate earlier 
for the modified evolution. 

Table H] compares the yield ratios from the resonance-decay model, within the region of 
optimal model-parameter values, with the available central Au+Au collision data from RHIC 
at y/s = 130A GeV The temperature of the resonance mass spectrum, Th, and the 
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TABLE I: Particle yield ratios from the resonance decay model compared to the RHIC Au+Au 
central-collision data at y/s = 130 A GeV. In the model, the starting baryon number of -Bo = 3 is 
assumed, and the mass of either too = 100 GeV or too = 95 GeV with a corresponding resonance 
mass spectrum temperature of either Th = 170 MeV or Th = 175 MeV is used. Results obtained 
when suppressing fusion reactions, for Th = 170 MeV, are provided in parenthesis. 



Yield-Ratio 


Decay-Model Results 


Experimental 


Collaboration Ref. 




T H = 170 MeV 


175 MeV Data 






p/7r- 


0.06 (0.05) 


0.07 


0.07 ± 0.01 


STAR 


[18] 


7T~ /tT + 


1.01 (1.01) 


1.01 


1.00 ± 0.02 


PHOBOS 


[19] 








0.95 ± 0.06 


BRAHMS 


[20] 


K+/K- 


1.073 (1.078) 


1.074 


1.092 ± 0.023 


STAR 


[21] 








1.28 ± 0.13 


PHENIX 


[221 








1.09 ± 0.09 


PHOBOS 


[19] 








1 12 ± 07 


RT3 AHMS 


[901 

IZUJ 


K~/ir~ 


17^ I'D 1 8^ 


0.175 


0.146 ± 0.024 


O 1 -Ti. 1 V 


I/O] 


p/p 


0.65 (0.63) 


0.66 


0.65 ± 0.07 


STAR 


[24J 








0.64 ± 0.07 


PHENIX 


f22] 








0.60 ± 0.07 


PHOBOS 


[19] 








0.64 ± 0.07 


BRAHMS 


[25] 


A/A 


0.72 (0.69) 


0.73 


0.71 ± 0.04 


STAR 


[21] 








0.75 ± 0.19 


PHENIX 


[26] 




0.81 (0.73) 


0.80 


0.83± 0.06 


STAR 


[2JJ 




0.82 (0.76) 


0.82 


0.95 ± 0.15 


STAR 


[2JJ 


A/h- 


0.054 (0.048) 


0.060 


0.054 ± 0.015 


STAR 


[27] 


TL/br 


0.040 (0.034) 


0.046 


0.040 ± 0.015 


STAR 


[27] 


A/p 


0.57 (0.58) 


0.59 


0.89 ± 0.22 


PHENIX 


[26] 


A/p 


0.63 (0.63) 


0.66 


0.95 ± 0.24 


PHENIX 


[26J 


Er/h- 


[6.0 (5.0)] xl0~ 3 


6.9x10" 


- 3 [7.9±l.l]xl0~ 3 


STAR 


[28] 


z + /h- 


[4.8 (3.8)] xl0~ 3 


5.6x10" 


- 3 [6.6±0.8]xl0~ 3 


STAR 


[28J 


S-/7T- 


[7.5 (6.7)] x 10-3 


8.9x10" 


- 3 [8.8±0.4]xl0~ 3 


STAR 


[29] 


S-/A 


0.108 (0.109) 


0.115 


0.193 ± 0.032 


STAR 


[28] 


I+/A 


0.118 (0.114) 


0.122 


0.219 ± 0.037 


STAR 


[28J 


n/br 


[2.7 (2.4)] xl0~ 3 


3.2x10" 


- 3 [2.2±0.6]xl0~ 3 


STAR 


[28] 
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initial baryon-number-to-mass ratio, Bq/ttlo, can be adjusted by simultaneously considering 
the ratios of p/tt + and p/p. An optimal agreement between the model and the central RHIC 
data is obtained for Tjj pa 170 MeV and B /m pa 0.030 GeV -1 . The model results tend to 
be only weakly sensitive to m : with an increase in mo a slight increase in Bq/uiq is favored, 
that can be traced to the factor f(m — m g ) in the density of states. 

The overall agreement between data and the decay-model calculations in Table U is quite 
remarkable, given that only two parameters, Tjj and Bo/mo, are adjusted. One should note 
that the optimal Hagedorn temperature of Tjj = 170 MeV is close to the critical temperature 
of T c ~ 170 MeV for a transition to the color-deconfined QGP phase obtained in the lattice 



QCD calculations at zero net baryon density [7(, and it is also similar to the chemical 
freeze-out temperature T^k = 174 MeV extracted from the analysis of RHIC data within a 
grand-canonical model 3j]. 

At the general level, the calculations are quite good in reproducing yield ratios involving 
strange particles. A more detailed examination, however, reveals some potential deficit of 
multistrange baryons and antibaryons. Possible reasons for the deficiency, to be investigated 
in the future, include: the possible role played b y th e non-resonant strangeness-exchange 
processes and by the multiparticle processes p\ and, further, the possible sensitivity 
of strangeness production to an early system dynamics 32( and, specifically, to strangeness 
fluctuations for early resonances and/or to details in flavor- dependence of state density. 

Besides the yield ratios in the standard model evolution, the corresponding ratios from 
the evolution with suppressed fusion processes are given in parenthesis in Table |J for Tjj = 
170 MeV. Though the suppression of fusion alters particle abundances early on in the system 
development, the final yields turn out to be rather similar, quite uniformly across the particle 
species. Only a careful examination reveals that the fusion suppression enhances slightly the 
production of pions and other light mesons and reduces slightly the production of heavier 
strange particles. 

Within a moderate range, there is no strong preference for a particular Hagedorn tem- 
perature. As Table H] shows, similar yield ratios are obtained for Tjj = 175 MeV as for 
Tjj = 170 MeV, if, in the context of the RHIC data, a slightly reduced initial mass of 
m = 95 GeV is assumed. 

In Fig. El we show the particle yield ratios that turn out to be particularly sensitive 
to the Hagedorn temperature under a given constraint. In one case, while varying the 
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FIG. 2: Dependence of particle yield ratios on Hagedorn temperature for a fireball characterized by 
the Bo/mo ratio adjusted to reproduce either the yield ratio of p/ir + = 0.11 (upper set of lines) or 
p/p = 0-65 (lower set of lines). The ratios chosen for the adjustment represent the central Au+Au 
data at s/s = 13CL4 GeV. 

temperature, we adjust the initial fireball's ratio Bq/trq to reproduce the ratio p/p ~ 0.65 
for the RHIC data (lower set of lines). In that case, the baryon-to-antibaryon ratios, B/B, 
remain rather stable with Tjj variation; the strongest variations are observed for the ratio 
of strange baryons to the negatively charged hadrons or to the negative mesons. On the 
other hand, if we adjust the ratio B /m to reproduce the ratio p/n + ~ 0.11 (upper set of 
lines), strong variations are observed for the B/B ratios. No matter what fitting strategy is 
followed, a reasonable agreement with the data is obtained within the Hagedorn temperature 
range of T H = 165 - 175 MeV. 

We next confront our resonance decay model with the SPS abundance data from the 
central Pb+Pb collisions at the laboratory energy of 158A GeV. As illustrated in Fig. 01 
an optimal agreement with the data is obtained for mo = 100 GeV when assuming (at 
Tjj = 170 MeV) a starting baryon number of B = 26. While the general agreement is 
rather good, we note that the calculated Q + /Q~ ratio is about 60% larger than the data. It is 
likely that the assumption of a larger number of lighter initial resonances would improve the 
agreement; in the thermal model the discrepancy is tauted as strangeness undersaturation j^j]. 

Measured kinematic spectra of particles from central heavy ion collisions exhibit the 
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□ Decay Model 



FIG. 3: Particle yield ratios from the resonance decay model compared to the SPS Pb+Pb central- 
collision data at the laboratory energy of 158A GeV. In the model, at Th = 170 MeV, the starting 
baryon number of Bq = 26 was assumed, in combination with the mass mo = 100 GeV and 
strangeness Sq = 0. 



effects of collective expansion. As may be expected, with suppressed fusion processes, our 
model produces kinematic spectra characterized by slope temperatures close to T H (see 
also Ref. In the standard version of the model, the sequences of decay and fusion 

generate some collective motion, but not enough to explain the transverse RHIC spectra. 
For pions e.g. the slope temperature raises by 4% compared to the version without fusion. 
Moderate increases in the elastic cross section, such as to an overall 10 mb, raise kinematic 
temperatures further, 6% for pions, but not well enough to approach data. As a consequence, 
it is necessary to assume the presence of some collective motion early on in system evolution, 
leading to resonances that exhibit space-momentum correlations. It might be that the 
dynamics, beyond statistics, needs to be involved in the predominant decays, involving the 
interior degrees of freedom of resonances. The first resonances might also emerge at finite 
transverse velocities. In either case, degrees of freedom beyond resonances would be involved. 

The transverse mass spectra displayed in Fig. 0] are obtained by folding a common collec- 
tive velocity field with the spectra from our decay model. Specifically, we assume a uniform 
transverse velocity distribution, d 2 N/d[3^ = 0(/3 t — (3 m3X ). The spectra for all the hadrons 
can be best described, at Th = 170 MeV, with a uniform velocity field of /9 max ~ 0.40, 
corresponding to an average flow velocity of {(3 t ) = 2/3 max /3 ~ 0.27. Notably, much less 
early flow, as characterized by {(3 t ) ~ 0.14, is required to explain the particle spectra at SPS 
energy, see Fig. El 
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FIG. 4: Transverse mass spectra of midrapi dity hadrons at RHIC. Solid symbols are the data 
fromPHENIX 3 

(left panel) and STAR |27l.l28j| (right panel) collaborations from central Au+Au 



collisions at the energy of y/s = 130A GeV. The lines represent the model spectrum from convo- 
luting the spectrum of the resonance decay model, at Tjj = 170 MeV, with a uniform transverse 
collective velocity distribution characterized by ((3t) = 0.27. 
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FIG. 5: Transverse mass spectra of midrapidity hadrons at SPS. Solid symbols are the data 



from NA44 3^] (left panel) and WA97 [36( collaborations from central Pb+Pb collisions at the 
energy of 158 AGeV. The lines represent the model spectrum from convoluting the spectrum of the 
resonance decay model, at Tjj = 170 MeV, with a uniform transverse collective velocity distribution 
characterized by (/3t) = 0.14. 



In summary, we have formulated a statistical model of hadron resonance formation and 
decay. Within the model, the density of hadronic states in mass is described in terms 
of a universal Hagedorn-type temperature. We have demonstrated that both the RHIC 
and SPS abundance data can be suitably described in terms of resonance decays at the 
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spectral temperature of T# ~ 170 MeV, even when pursuing the extreme assumption of a 
single heavy resonance populating the investigated rapidity region. To explain the data for 
particle spectra, we needed to invoke additional collective motion beyond that generated in 
the hadronic interactions. 
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